Tunneling through Al/AlOx/Al junction: analytical models and first principles simulations 
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We study from first principles the transport properties of Al/AlO^/Al tunnel junctions. On this basis, we 
analyze the reliability of two analytical models for the conductance, namely the trapezoid potential barrier model 
and a tight-binding model. Our findings show that (;') the interface width used in the models is determined by 
the electronic density profile, and it is shorter than the width one expects from the atomic arrangements; (ii) the 
effective mass, found to be about on third of the free electron mass, can be determined from the oxide band- 
structure calculations, and (Hi) the barrier height is given by one fourth of the bandgap in the oxide, which 
explains the apparently small values found for these junctions experimentally. 

PACS numbers: 73.40.Rw; 73.40.Gk; 71.15.-m 



I. INTRODUCTION 

Tunneling of electrons through alurninum-aluminum oxide 
(Al/A10 r /Al) junctions is one of the prototypical examples of 
quantum-mechanical tunneling in solid state physics. Fisher 
and Giaever in their pioneering work 1 demonstrated the tun- 
neling character of the transport of electrons through this in- 
terface and by comparing their results with the predictions of 
Holm 2 for tunneling through vacuum gap, they initiated the 
interpretation of tunneling measurements through thin metal- 
insulator-metal junctions using potential barrier model. The 
minimal form of this model contains two parameters - the bar- 
rier width d, which indicates the physical width of the oxide, 
and its height W, given by the energy difference between the 
Fermi energy and the bottom of the conduction band in the ox- 
ide. In practice, several other parameters enter the model 3 -" 5 ^ 
the electron's effective mass in the oxide or the dielectric con- 
stant of the oxide used within an additional image-charge po- 
tential. Further parameters are used for fine-tuning the shape 
of the barrier, e.g. it's asymmetry 5 . Clearly, having a large set 
of parameters, it is no surprise that the simple barrier model 
can be fitted to the experimental current-voltage characteris- 
tics welP^l, but at the same t ime, it rises questions about the 
relevance of the model itself LL2UU For example, the inclusion 
of the image potential can have a significant effect on the ef- 
fective barrier width, but its presence depends on the time 
scales of the tunneling electrons and the interface plasmons 
in the metafSEH. 

On the other hand, much more detailed and parameter- 
free models of the interface can be constructed using first 
principles calculations^^ even though the size of the mod- 
eled interfaces is somewhat restricted due to the numerical 
cost of these ca lculations. Nevertheless, in many experi- 
men j JHIlinnillQlE] jjjg sni( jj ec [ interfaces have widths within 
the reach of ab initio simulations so that the accuracy of the 
potential barrier model to the interpretation of tunneling data 
can be tested. Specifically, Jung et af^lpresented such a study 
comparing the character of the equilibrium projected density 
of states of the Al/A10 r /Al interface obtained by a first princi- 
ples simulation with the potential barrier model. They found 



that the parameters of the potential barrier model fitted to the 
experimental data are in qualitative agreement with the param- 
eters of the first principles calculations. The potential barrier 
model included the image potential and hence also the dielec- 
tric constant which effectively narrowed and lowered the po- 
tential barrier. 

In this work we test the performance of the potential barrier 
model by comparing the predicted conductance to ab initio 
calculations^ 2 ^. We test this for A1/A10 A /A1 junctions of 
four different widths d and show that it is essential to use an 
effective mass in the oxide and an effectively shorter width of 
the tunneling region within the potential barrier model. We 
also present an analytical tight-binding model for the conduc- 
tance that describes the ab initio results more accurately than 
the potential barrier model. The parameters of the latter are 
extracted from the ground state ab initio calculations of the 
junction. In Sec.|Il]and|nI]we introduce the analytical details 
of the models. The ab initio results for ground state proper- 
ties of the studied junctions are presented in Sec. IV and [V] 
together with the computational parameters used in the cal- 
culations. Finally in Sec. |Vl]we compare the conductances 
obtained using the ab initio calculations and the conductances 
obained from the analyical models. 



II. POTENTIAL BARRIER MODELS OF THE INTERFACE 

The starting assumption of the potential barrier model is 
that inside the metallic electrodes, on the left and right of the 
insulator, the electrons behave like free quasi-particles with 
their energy being in a separable forrrP^ 



E=E-+En = k 2 v /2 + kl/2 



(1) 



where k z is the component of electron's momentum perpen- 
dicular to the interface and fcp the component of momen- 
tum parallel to the interface. The current density, induced 
by an infinitesimal bias voltage, consists of a sum of con- 
tributions from the electrons occupying states in the energy 
window around the Fermi energy Ep, with their momentum 
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FIG. 1: (Color online) The trapezoid potential barrier used to model 
ultra-thin AlOx interface. 



opposite to the drop of the bias voltage (k z > 0). Hence, the 
conductance per area is given by the expression^! 
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where T(E Z ) is the transmission function, the probability for 
an electron to pass through the junction, and Ep is the Fermi 
energy. 

The simplest expression for the transmission T(E Z ) is based 
on a metal-vacuum-metal interface^! where the barrier 
height W is given by the potential energy in the vacuum with 
respect to the Fermi energy of the metal. 

In Sec.[VT]we will demonstrate that there are two essential 
features of the potential barrier models that need to be taken 
into account for the description of ultra-thin interfaces: (1) the 
barrier needs to have transition regions between the metal and 
the insulator of width Ad, where the potential energy changes 
continuously, (2) the effective mass of the electrons in the in- 
sulator needs to be accounted for. These two requirements 
can be fulfilled by using a specific shape of the potential bar- 
rier. In this work we use a trapezoid potential barrier (TB) as 
defined in Fig.[T] 

On the other hand, in Sec.[VT]we will also demonstrate that 
approximate expressions for the transmission as well as for 
the energy integration in Eq. Q are sufficient for an accurate 
evaluation of the model conductance. For the trapezoid po- 
tential barrier , the WKB approximation for the transmission 
gives 



T(E-E\\) = esxp{-F(£-JS||)} 
F(E — En) = 



2m eS [W(z)+Eu]dz, 



(4) 
(5) 



where m e ff is the effective mass of the electrons in the insula- 
tor, — d{E\\) and d(E\\) give the region where [W(z] + E\\] > 0, 
and W(z) is the trapezoid potential profile. Accounting only 
for the largest contribution from the states close to the Fermi 
energy in the integral in Eq. Q, E^ ~ Ep, we obtain the 



following simple anaytical expression: 
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(7) 



F'(E F ) = - m (d w +2M EF ). (8) 

We will refer to Eq. ^ as the TB A ™ eff model (A stands for 
"analytical" as compared to the numerically calculated trans- 
mission for the trapezoid potential barrier - TB^). We note 
that the introducion of two transition regions of width Ad adds 
to the exponent of the transmission amplitude only a small 
fraction of Ad, namely (2/3)AdE F - This results in a sub- 
stantial increase of the conductance which is needed for the 



agreement of the TB model and ab initio results (see Sec. VI I. 



III. ATOMIC sp MODEL OF THE INSULATOR 

It is typically assumed that the barrier height in the poten- 
tial barrier model corresponds to energy distance between the 
Fermi energy and the closest among the valence or conduc- 
tion bands of the insulator, or even to its whole bandgap. How- 
ever, fits of the potential barrier model to experimental data 
often lead to unphysically small values if one follows this in- 
terpretation. Various arguments like interface roughness^ or 
image potential 4 have been suggested to correct for this un- 
derestimation, but perhaps the most important one - the prin- 
cipal difference in the energetic spectrum of the real insulator 
and the vacuum gap - received less attentiorPSO 

To account for a more realistic electronic structure of the in- 
sulator we consider a minimal tight-binding model of a sp-like 
insulator with rock-salt crystal structure. For our purposes, 
the cation with s-like orbital plays the role of aluminum and 
the anion with /?-like orbital the oxygen atom. While this 
is different from the true structure of alumina, this model 
works surprisingly well even for the disordered aluminum ox- 
ide found in our interfaces, as will be shown in Sec.|IV| 

The sp model has four parameters: the onsite atomic en- 
ergies of the cation (£ s ) and anion {e p ), the hopping matrix 
element between the two atoms (f), and the length of the edge 
of the conventional unit cell (cube) a. A standard calculation 
leads to valence (v) and conduction (c) band energies 



£ c /v(k) =Ef ± 
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(9) 



where m e ff = E s /(2t 2 a 2 ) is the effective mass of the elec- 
trons close to the conduction band minimum, equal in magni- 
tude that of the valence band maximum. The two bands are 
separated by the bandgap E g = e p — E s , and the energy in the 
middle of the gap is 



(10) 



In the tunneling regime, the current is carried by the elec- 
tronic states in the bandgapj 27 " 29 ! i.e. the evanescent Bloch 
states with imaginary wavenumber k z = iK\ 

^ || (r)~ e -^ e ik ir r MK ,k || (r). (11) 
The WKB-like result for the transmission takes then the form 

T°P(E)~\<j> K ^(d)\ z ~e- 2 « E ' k \O d , (12) 

where d is a vector normal to the interface with the length 
given by the width of the interface (|d| ~ d). K(E,k« ) can be 
obtained from Eq. (|9]l using the substitution k z = iK therein. 
The transmission can be then used in the calculation of the 
conductance in Eq. [2] The largest contributions to the con- 
ductance come only from Ka/2 < 1, kna/2 < 1, so that the 
sin( ) functions in the dispersion can be expanded in Taylor 
series. Keeping only the first two terms we fincP 




FIG. 2: (Color online) The 4L structure (above) and the correspond- 
ing averaged electronic density of the occupied transmitting states. 
The dashed lines indicate the positions of the metal-oxide boundary 
obtained according to Eq. (TTSl. 



K(E,h) 



v{E)m eS E g /2 + kl + l&, (13) 



2[v(E)m eS E g /4 + E ]l \, (14) 



where we have introduced a multiplicative factor v(E) ac- 
counting for the relative distance of the energy E from the 
middle of the gap, 



v{E) = 1-4 



E-El 



(15) 



which is close to 1 for E ~ £>. We note that by using the 
Taylor expansion the model becames independent of the size 
of the conventional cell a. The transmission T k p (E) is simi- 
lar to the WKB result for a potential barrier [Eqs. (|4]|5]) for a 
constant barrier height W]. Hence, making the same approx- 
imations as in Sec. [n] and substituting W — > v(Ef)E g /4 we 
find an analytical expression for the transmission through a 
sp insulator of width d precisely of the form of Eq. ([6j, where 



F sp {E F ) = 2^v{E F )m ef[ E g /2d, (16) 

F '^ Ef) = A,f F f W h d - (17) 
y/v(E F )m ee Eg/2 

This represents one of the main results of our paper: the po- 
tential barrier height W is related to the bandgap through the 
relation W = v(E F )E g /4, Since the Fermi energy in our junc- 
tions is close to the center of the gap (Sec. IV i where we 
have v(Ep) ~ 1, we expect that the bandgap is about four 
times larger than the barrier height obtained from the fits to 
the experimental data. This explains the typical situation in 
Al/AlOj/Al junctions where W can be as small as 2eV or less, 
which is to be compared with the bandgap of alumina be- 
ing about 7 — 9eV. Further comparisons will be made in the 



Sec. VI where the sp model is compared to the ab initio cal- 
culation of the conductance. 



IV. FIRST PRINCIPLES CALCULATIONS OF THE 
AL/ALO A /AL INTERFACES 

The Al/A10 x thin film is well k nown for its difficulties to 
be grown in an ordered 

formSED. 

The process of oxidation 
consists of a quick chemisorption of oxygen on a clean Al sur- 
face which is followed by a complex diffusion process lead- 
ing to vari ous widths of the interface which is typically disor- 
dered™!!] 

The model that we consider is on the other hand 
relatively simple and ordered. We followed Jennison 35 ! 36 ! at 
constructing chemisorbed layer of oxygen on an ideal Al(l 11) 
v3 x v3 surface (three Al atoms per layer), modelled as a 
slab 6 layers thick (left electrode). Next we were adding Al 
and O atoms and relaxed the geometry until we found a stable 
interface having two layers of oxygen atoms (2L). Finally we 
enclosed the interface with four ideal Al(lll) layers (right 
electrode) and connected it with the left electrode through pe- 
riodic boundary conditions. Performing this procedure two 
different geometries of the interfaces were identified: (1) an 
asymmetric structure, corresponding the the ultra-thin A10 A 
layer investigated by Jennison, and (2) an symmetric struc- 
ture which did not contain the layer of chemisorbed oxygen 
next to the bottom Al electrode. More details on the differ- 
ences between the asymmetric and symmetric structures can 
be found elsewhere 3 -^, in our present work we will consider 
only structures derived from the asymmetric geometry. 

Motivated by the geometry of the asymmetric 2L interface 
model we have constructed thicker A1/A10 X /A1 by adding one 
(3L), two (4L) or three (5L) full oxygen layers sandwiched 
between monoatomic (All) or diatomic (A12) layers of alu- 
minum. The resulting geometries were optimized until the 
forces on the atoms were smaller than 0.002 Ha/ag, while 
the Al atoms beyond the first layer of bulk metal were kept 
fixed. An example of the resulting geometric structure of 4L is 
shown in Fig. [2] We should mention that these models are not 
necessarily the only ones possible for the interface of the con- 
cerned width. Due to the above described tendency of A10 x 
systems towards disorder, we expect that many different vari- 
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System 


d[A] 


Eg [eV] 


AEf [eV] 


v{E F ) 


2L 


4.5 


7.0 


1.5 


0.82 


3L 


5.5 


6.5 


-0.25 


0.99 


4L 


7.8 


6.5 


-1.0 


0.91 


5L 


9.8 


6.5 


-1.0 


0.91 



TABLE I: Values of the interface widths, band-gaps, Fermi energy 
shifts and the shift factor v(Ep) obtained from ab initio calculations. 



ations could be found with larger surface cells. The structures 
identified here need to be taken as few samples of the great 
variety of possible geometrical arrangements. However, the 
comparison of the projected density of states (PDOS) for sym- 
metric and asymmetric 2L junctions (see Ref. 11371 ) suggests 
that these differences lead to small changes in their conduc- 
tances. 

All of the ground state properties and optimizations were 
done using the Quantum Espresso distribution 38 . We have 
employed the PBE exchange-correlation functional, atomic 
cores were described using ultra-soft pseudopotentials result- 
ing in well converged electronic structure close to the Fermi 
energy, using a cutoff energy of 12.5 (125) Ha for wavefunc- 
tions (charge density). Due to the large size of the supercell, a 
6x6x1 Monkhorst-Pack k-point grid was sufficient to con- 
verge the total energy and the electronic density. 

In Fig. [2} in parallel with the geometric structure of the 4L 
interface, we show the profile of the plane-averaged electronic 
density of the scattering states An(z) (e.g. localized states on 
oxygen atoms are not included). We see that the rapid drop 
and increase in the density appears at the boundary between 
the metal and the oxide. We use An(z) as the quantity for the 
determination of the interface width from our ab initio calcu- 
lations, in close analogy with the determination of the position 
of surfaces at metal-vacuum interfaces 39 ; for the left bound- 
ary we use 
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(18) 



where zi is a position in the center of the insulator. Similar 
expression is used for the determination of the right boundary 
Zr which together with zl gi ve the estimate of the interface 
width d = zr— Zl used within our potential barrier and sp 



models in Sec. VI The resulting interface widths are given in 
the Table [j] In the following we will also refer to the width of 
the transition region beween the metal and the oxide, which 
can be estimated from the averaged density to be An « 2.0A. 
This value will be used for the determination of the width of 
the transition region in the potential barrier model (Fig. |TJ. 

We note that for the calculation of the positions zl/r we 
do not necessarily need to use the density An(z) obtained 
from the scattering states in the transport calculations (see 



Sec VI I, but is it equally good to use the partial density of the 
states close to the Fermi energy that can be obtained from any 
ground state code (e.g. Quantum Espresso). On the other 
hand, the total electronic density or the Kohn-Sham potential 
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FIG. 3: (Color online) The comparison of the average Kohn-Sham 
potential v^(x), the total density n(z) and the contribution to the 
density from the states close to the Fermi energy An(z). Clearly, the 
latter can be unambiguously used for the definition of the width of 
the interface using the Eq.|18| 



(which is frequently but incorrectly believed to be the origin 
of the potential barrier in the model from Sec. [Ilji are not suit- 
able for this calculation, as it is clearly demonstrated in Fig. [3] 



The second important parameter of the potential barrier and 
the atomic sp models is the insulator band gap E„. It can be 
extracted from the projected density of states (PDOS), where 
the Kohn-Sham eigenstates of the interface are projected on 
atomic orbitals. Figure [4] shows the PDOS for the 4L inter- 
face, where the PDOS of atoms in each layer are added to- 
gether, green lines corresponding to the Al layers and the red 
lines to the oxygen layers. The oxide bandgap is estimated 
as the energy distance between the onset of the valence bands 
on the oxygen atoms below the Fermi energy, and the onset 
of the mixed Al and O bands above the Fermi energy. From 
the PDOS we can also obtain the energy distance between the 
midgap energy and the Fermi energy, AEf, needed for the sp 
model. The calculated bandgaps and AEp for all studied 
interfaces are collected in the Table [I] Interestingly, in spite 
of the well known bandgap problem of the DFT 40 " 42 , these 
bandgaps appear to be in very good agreement with recent 
experimental r esults for the AI/AI2O3 interfaces^ 4 ! which 
found E g = 6.4eV. 

While the bandgap stays roughly the same for all of the 
studied interfaces 2L — 5L, the Fermi energy shifts with re- 
spect to the middle of the gap from positive (conventionally 
called the electron tunneling regime) to negative values (hole 
tunneling). However, the factor v(Ep) stays close to one in 
all the cases (see Table|I]), as anticipated already in Sec|ll] The 
energy difference between the bottom of the conduction band 
and the Fermi energy determined experimentally 4 - 3 -' was found 
to be E c — Ep = 2.9 ± 0.2eV which is leV smaller than the 
DFT value found here for 4L and 5L, but on the other hand, 
in good agreement with 2L and 3L, which perhaps indicates 
larger sensitivity of this quantity on the particular system. 





-10 -8 -6 -4 -2 2 4 6 8 10 

E - E F [eV] 

FIG. 4: (Color online) The projected density of states for the AlOx 
4L interface. The distance of the valence band to the Fermi energy, 
being less than half the distance to the conduction band, is taken as 
the effective barrier height. 
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FIG. 5: (Color online) The imaginary (left) and real (right) band 
structures from ab inito calculations compared with the band- 
structures of the sp and potential barrier models. The sp model 
gives an excellent fit for both real and imaginary band-structure for 
m eff = 0.35 and E g = 4.0eV. 



V. ELECTRONIC STRUCTURE OF AN IDEAL 
INSULATOR 

The potential barrier model as well as the sp model also 
rely on the knowledge of the effective mass m e g of the elec- 
trons in the insulator or barrier region. To calculate it we have 
considered a first principles model of the insulator extracted 
from the geometry of the 4L junction. Namely, it consists 
of a supercell of length I = 8. 1 las in the z direction and with 
identical dimensions in the two remaining in-plane directions 
[i.e. \/3 x \/3 Al(lll)], The supercell contains two layers of 
oxygen and two layers of 2/3 filled Al planes, (the 3rd and 
4th oxygen layers in Fig. [2] from the left and their immedi- 
ately following Al layers respectively). This way, the chemi- 
cal composition actually corresponds to alumina, AI2O3. 

The DFT ground state calculation has been done with the 



same specifications as for the full interface (Sec. IV 1 except 
for the k-point grid being here 6x6x6 due to smaller extent 
in the z direction. The following band-structure calculation 
has been done using the PWCOND prograrrJS! that is capable 
of obtaining the so called complex band-structure, i.e. en- 
ergy bands for imaginary as well as real Bloch fc-vectors. We 
have checked that calculations of the band-structure for real k- 
vectors using the Quantum Espresso and the PWCOND gave 
identical results so that the parameters involved in the PWCOND 
program were correctly chosen. 

The band-structure along the direction normal to the in- 
terface (z) is shown in Fig. [5] First of all we note that the 
bandgap obtained here, E™ ~ 4eV (in agreement with the pre- 
vious DFT-PBE results for bulk y-Al 2 o|2E0), is significantly 
smaller that the bandgap extracted from the PDOS of the full 
junction (~ 6.5eV). Interestingly, the experimental value of 
this phase of alumina is E e g xp = 7eV, which can be obtained 



also computationally if the DFT-PBE result is followed by a 
GW calculatiorP. 

The DFT band structure in Fig. [5] is fitted with two model 
dispersions. The TB model uses a free-electron like disper- 
sion E c (k) = e c +k 2 / (2m e ff) which after fitting gives the ef- 
fective mass m e ff = 0.35. The atomic sp model [Eq. |9])] in 
the approximation ka/2 < 1, which is used in the analytic ex- 
pression for the conductance, gives (for k x = k y = 0) the dis- 
persion: 



E c/v (k) 




(19) 



The parameters of the fit given in Fig. [5] are m e ff = 0.35, 
Eg = 4.0eV and Ep = leV. Our value of the effective mass is 
to be compared with the electron's mass obtained from DFT 
calculations for ideal a-A^Cb crystal, m e ff « .4^, and fits 
to experimental / — V characteristics, m e ff ~ 0.23^H 

We see that the sp model works very well for real as well 
as imaginary band-structure close to k = 0. While both mod- 
els give the same effective mass, the values of K for the free- 
electron like dispersion are larger by ~ 50% (as indicated by 
arrows) which contributes to prediction of smaller conduc- 
tances within the potential barrier model given the interface 
width is the same, as will be shown in the following section. 



VI. THE CONDUCTANCE 

Transport properties of the junctions were obtained us- 
ing the transfer matrix method 49 implemented in the PWCOND 
code^J , using plane-wave basis and ultra-soft pseudopoten- 
tials. For the given self-consistent Kohn-Sham potential (ob- 
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code 


2L 


3L 


4L 


5L 


PWCOND 


0.109 


0.0166 


0.00245 


0.000279 


WanT 


0.0668 


0.00730 


0.00224 


N/A 



TABLE II: Values of the conductances in multiplies of e 1 /h x A z , 
where A ; = 74. 23^ is the area of the supercell perpendicular to the z 
direction, calculated by the PWCOND and WanT codes. The differences 
are similar to the differences between the PWCOND results and the sp 
model. 
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FIG. 6: (Color online) The ab initio and model conductances. The 
square barrier (green dashed line), the trapezoid barrier with transi- 
tion regions AD / (pink dotted line) and the inclusion of the effec- 
tive mass (orange dash-dotted line) form a sequence of improvements 
of the potential barrier to model towards the true atomic system. Fi- 
nally, the analytic sp model (black full line) with parameters taken 
from ground state ab initio simulations gives very good agreement 
with the full first principles calculation of the conductance. 



tained from the ground state calculations, see Sec. IV i, the 
conductance was converged with respect to the fcp grid; go- 
ing from the 6x6 mesh (used for the presented results) to 
a 10 x 10 mesh the change in the conductance has been 
found to be < 5%. Furthermore, for testing purposes, the 
conductances for the 2L, 3 L and AL interfaces were also cal- 
culated using the WanT co de 24 25 , where a completely differ- 
ent method based on maximally localized Wannier functions 
is implemented. Results are reported in Table [H] and compare 
well with the previous set, though slightly underestimating the 
absolute values. 

In Figure [6] we show the dependence of the conductances 



per unit area on the interface width d, determined in Sec. IV 
in comparison with the two models considered in Sections]!!] 



and III The horizontal error bars accompanying the ab initio 
conductances, Ad s» 2A, indicate the width of the transition 
region between the metal and the insulator, which is taken 
from the averaged density profile, Fig. [2] 

First we consider the potential barrier model with effective 
mass equal to one, where the calculation of the transmission 
as well as its energy integration (Eq. [5]) are done numerically 
exactly (TB^). The potential barrier is of the form given in 
Fig. [T] where d w — d — Ad. The conductances are shown 



as the blue crosses, where the height of the energy barrier 
W = 0.5neV, n = 1,2,3,4,5 is increasing from top to bottom. 
The pink-dotted line is the conductance of the same potential 
barrier of width W = 2eV, but evaluated using the approximate 
formula [Eq. d6)]. As anticipated in Section[II] we see that in 
view of the overall differences, the approximate but analytic 
formula is very satisfactory and the numerical calculation of 
the transmission of its energy integration is not really needed. 

We see that in principle, we can achieve agreement between 
this model and the ab initio results if we choose W ~ 0.5eV, 
but this is in stark contrast with the estimates of the potential 
barrier height from the PDOS, typically taken as the distance 
between the Fermi energy and the nearest band in the oxide 
(e.g. the valence band in the oxide in 4L structure according 
to Fig. HJ, here expected to be W ~ 2eV. 

The green-dashed line is a conductance corresponding to 
a simple square potential barrier with W = 2eV and effec- 
tive mass equal to one, and we see that plain square barrier 
model goes in the wrong direction. The use transition regions 
of width Ad does shift the potential barrier model in the right 
direction, particularly for very short interfaces, where the ef- 
fective mass within the insulator does not seem to play an im- 
portant role. Hence, use of the transition region between the 
metal and the insulator of width Ad, given by the spatial extent 
of the drop if the electronic density between the metal and the 
oxide, is essential for the TB model. 

The red dash-dotted line gives the conductance accord- 
ing to Eq. (|6| with the ab initio determined effective mass 
m e ff = 0.35 and W = 2eV. The effective mass significantly im- 
proves the agreement of the potential barrier model with the 
ab initio conductance, while keeping the barrier at the "rea- 
sonable" value, motivated by offset between the Fermi energy 
and the valence band maximum. 

Finally, the full black line corresponds to the atomic sp 
model with the effective mass m e ff = 0.35, band gap E g = 
6.5eV and the barrier width dw = d — Ad. The use of this 
reduced width dw is motivated by two observations: (1) 
in Sec. [TT] we have seen that the linearly increasing poten- 
tial at distance Ad contributes to the exponent of the con- 
ductance [Eqs. d6fl7])] through a much smaller contribution 



Ad Ep = W/{W+E F )Ad ~ 0.15M. (2) in the TB model we 
have seen that the use of a shorter barrier, effectively given 
by dw +2/3 Ad [Eq. Q], is important to compare well with 
the ab initio conductances. Hence we expect that also in the 
sp model, the oxide width (i.e. the equivalent of the potential 
barrier) needs to be reduced almost to d — Ad, which is the 
value we use. As a result, the sp model is essentially on top of 
the ab initio conductances. While the improvement with re- 
spect to the potential barrier model with transition region and 
the effective mass is not that large, it is important to stress that 
the parameters of the sp model (E g , m e ff, d — Ad) correspond 
to the characteristics of the true ab initio model. 

It is interesting to attempt a quantitative comparison be- 
tween experimentally determined barrier widths and heights, 
and our ab initio and sp model results (Fig. [7]). As mentioned 
already in the introduction, there are experimental junctions 
that are now accessible to first principles simulations. Based 
on the rather unsatisfactory state of affairs in Fig. [7] we sus- 
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FIG. 7: (Color online) The ab initio conductances compared to se- 
lected experimental results, (a) Jung22 (b) GlooJ^l (c) HolmqvisP2l 
(d) Brinkmar. 5 , and recent experiments by Schaefer 9 . The model 
gives fairly rigid prediction of the conductances, even using the 
bandgap of the 0!-Al2O3, E g ~ 9eV. The likely source of these dis- 
crepancies is the experimental determination of the interface width. 



ultra-thin A1/A10 A /A1 junctions. We have compared atomistic 
first-principles calculations using the DFT-PBE framework 
combined with the Landauer formula, with the conductances 
obtained from a potential barrier and a tight-binding sp an- 
alytical models. We have shown that the expression for the 
conductance of the atomic sp model has the same form as that 
from the potential barrier model if the barrier height W is ex- 
changed for v(Ep)E g /A with v(Ep) ~ 1, which explains the 
small values of W obtained frequently in the past by fitting the 
potential barrier model to the experimental I — V curves. The 
accuracy of the analytical models has been tested by using 
parameters derived from ground-state DFT calculations. We 
have found that the oxide is characterized by effective mass 
tn e s = 0.35 and bandgap E g = 6.5eV. When these parameters 
are used in combination with the sp model, excellent agree- 
ment with the numerically calculated conductances is found. 
The interface width used in the models has been shown to 
correspond to the width of the well-developed oxide which is 
shorter by about Ad ps 2 A compared to the geometric width of 
the interface d. 
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